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SECTION I 
INTRODUCTION 


The AZUSA Geometric Dilution of Precision (GDOPS) Program (AZGD2) is 
designed to compute GDOPS for a standard KSC theoretical trajectory based on an 
AZUSA (or AZUSA-like) system. Ir, ternary the AZGD2 program uses a right-handed 
cartesian coordinate system located at the origin (%, M wit * 1 * downrange at 
some azimuth from north, Z crossrange 90 degrees greater than the azimuth of X, and 
Y perpendicular to the XZ plane. The output may be in this coordinate system or in 
the so called Apollo Saturn Coordinate System, which differs from the afore- 
mentioned system only in that its Z is downrange, its Y is crossrange, and its X is 
perpendicular, i.e., the axes are relabeled. 


Symbol 

a 

b 

AZ 

2 

e 

. <M,h 



a 


in 


/ 


°m 


X/Y ,Z 
X,Y ,Z 


a 


X 




Y 


CT z, a z 


u v 

*i,vj 


SECTION II 

DEFINITION OF SYMBOLS AND TERMS 

I 

Definition 

i 

Semimajor axis of Earth ellipsoid (input) 

Semiminor axis of Earth ellipsoid (input) 

Azimuth of X-axis (input) 

Square of ellipsoid eccentricity 

Geographic location (latitude, longitude, altitude) subscripts: 
o : origin, J : baselines, J = 1,5 (see below for identification) 
(input) 

Sigma for L-cosine, sigma for L~cosine rate (input) 

Sigma for M -cosine, sigma for M -cosine rate (input) 

Sigma for range, sigma for range rate (input) 

Vehicle cartesian position (input) 

Vehicle cartesian velocity (input) 

Expected X-position and X-velocity errors 

Expected Y-position and Y -velocity errors 
Expected Z-position and Z-velocity errors 

j 

Expected vector resultant velocity error 
Cartesian location of baseline points 


Symbol 


Definition 



L baseline end points 


M baseline end points 


5 


Range tracking point 



SECTION III 

COMPUTER PROGRAM EQUATIONS AND LOGIC 


A. GEOXYZ 


Geographic to cartesian position subroutine: 


N = 
o 


' 2 . 2 ^ 
1 - e sin cp 


U = (l\i “Hi ) cos r -P 
0 0 0 0 


V = [(1 - e ; N + h ] sin c p 
o L o o J o 


W = 0 
0 




, 2 . 2 
i - e sin cp. 


U , — (N , + h ,) cos C P, cos (X - A.,) 
J J J J o J 


Vj = [(i - e ; Nj + hj] sin <P 


W. = (N , T h .) cos f P.sin(A - X ) 
J J J J o J 


U = U , - U 
J o 


V = V . - V 
J o 


w = w. - w 

J o 
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X J “ “ u 

cos AZ cos^o *PV cos AZ sin*? 

Yj = U cos f P 0 + V sin Vo 

Zj = U sin AZ sin Vq - V cos sin AZ 

MAIN PROGRAM 

XLL = 

X 1 +X 2 
2 

YLL = 

Y 1 ' HY 2 
2 

ZLL = 

Z l +Z 2 

2 

XMM = 

x 3 +x 4 

2, 

YMM = 

Y 3 + Y 4 
2 

ZMM = 

Z3 + Z^ , 

2 

XLX = 

x ~ X 
a 2 1 

XLY = 

Y - Y 

Y 2 Y 1 

XLZ = 

Z - Z 
^2 1 

XMX = 

*4 " X 3 

XMY = 

Y 4~ Y 3 


~W Sin AZ 


-W cos AZ 


XMZ = Z 4 - Z 3 

B l = (XLX 2 +YLY 2 +ZLZ 2 ) 1/2 
B.. = (XMX 2 H-XMY 2 -i-XMZ 2 ) 1 ^ 

IVI 


CXL 

— 

XLX/B. 

CYL 

— 

YLY/B, 

CZL 

— 

ZLZ/B L 

CXM 

— 

XMX/B m 

CYM 


YMY/B m 

CZM 

— 

ZMZ/B m 

XXL 

== 

X - XLL 

YYL 

— 

Y - YLL 

ZZL 

= 

Z - ZLL 

XX M 

— 

X -XMM 

YYM 

— 

Y - YMW! 


. ZZM = Z - ZMM 


YYR = Y - X K 
5 

ZZR = Z -X, 

5 

DLMXYZ 

Partial derivative computation subroutine: 

(A = XXL, XXM; B = YYL, YYM; C = ZZL,ZZM; D ~ CXL, CXM; 

E = CYL , CYM; F~CZL,CZM; input arguments) 

wP wlP are elements of a row vector. (See PH1GAM, paragraph E.) 

X / j J 


p 

,.2 , d 2 . .2.1/2 
(A + B rC ) 

cvx = 

A/P 

CVY = 

B/P 

CVZ = 

C/P 

< 

(1 - CVX 2 ) D/P - CVX -CVY- E/P - CVX-CVZ-F/P 

- 

-CVX ’CVY' D/P + (1 - CVY 2 ) E/P -CVY -CVZ- F/P 

w™ = 

-CVX -CVZ- D/P - CVY- CVZ- E/P + (1 - CVZ 2 )-F/P 


DRXYZ 


Partial derivative compulation subroutine; 


R 


exxir i YYR‘ 


+ ZZR 2 ) 1/2 


W< R) = XXR/R 
W 2 R) = YYR/R 


Wg R) - ZZR/R 

<W< R) wiP Wf are elements of a row vector,) 


PH I GAM 


Partial derivative matrix subroutine: 



P 11 

p 12 

p 13 _ 


W i U 

W ( 2 L) 

W-> 

3 

p = 

P 21 

CM 

CM 

CL 

p 23 

— ' 

w m 
1 

W (M) 

W <M) 


P 31 

— 

P 32 

p 33 


W (R) 

1 

W (R) 

2 

wf 


VARCOV 


Variance-covariance matrix, and sigma computation subroutine: 

(S^ Sp = a ni/ = cr^ ; input arguments . Output 

i t 

arguments SX = Oy rs^} SY ~Oy Oyi SZ = o^) 


AL 

p n /s i ; 

BL - 

p i2 /s r 

CL 

?. n A* 

1 13 r "1 

AM 

P 21 /S 2 ; 

BM - 

P 22 /S 2 ; 

CM 

'■ P 23 /S 2 

AR 

■? P 31 /S 3 j 

BR 

P 32 /S 3 ; 

OR * 



AA = AL 2 + AM 2 + AR 2 

AB = AL'DL + AM • BM + AR • BR 

AC = AL-CL + AM-CM + AR-CR 

BB = BL 2 + BM 2 + BR 2 

BC, = BL.-CL + BM'CM + BR'CR 

CC = CL 2 +CM 2 +CR 2 


— 

h n 

h 12 

h 13 


AA 

AB 

AC 

h 21 

h 22 

h 23 

= 

AB 

BB 

BC 

h 31 

CM 

CA 

h 33_ 


AC 

BC 

CC 


H (1) = H" 1 , SX / SY =\/h^ , SZ 

(Matrix inversion accomplished by MINVRT subroutine, paragraph H.) 


G. S1GVEL 


Vector velocity sigma subroutine: 
V = (X 2 + Y 2 + Z 2 ) 1/2 


*2 ,( 1 ) 
= X 'll 
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DA 


DAB 


H. 


2X-Y’h 


( 1 ) 


12 


DAC 


2X • Z • 


DB 


* 2 ,( 1 ) 

Y h 22 


DBG = 2Y • Z • h 


( 1 ) 


23 


DC 


;. 2 ,( 1 ) 

z h 33 


= (DA + DAB + DAC + DB + DBC + DC) 1/2 /V 


MINVRT 

Matrix inversion subroutine: 




h n 

h 12 

h 13 

* 

H 

— 

h 21 

h 22 

h 23 




h 31 

h 32 

h 33 


A 1 

— 

h ll 

(h 22 h 33 “ h 23 h 32 

CM 

< 

— 

h 12 

(h 2 i h 

- h h 
33 n 23 31 

A 3 

= 

h 13 


32 “ h 22 h 31 


df 


H 


= a 1 -a 2 +a 3 

- 0, matrix is singular and inversion is impossible.) 
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* »i'*i juuv.wi y-y mrATy'i Wff nti • ar **?? ■"■.'If' ‘r"" v» ,,f i *• ** w t * 


H 


(C) .. 


(C) 

'll 


A l /h ll 


(C) - 
12 " 


(h 12 h 33 


h 32 h 13 5 


h 13 ~ h 12 h 23 ~ h °° h 


22 13 


l/C) = ~A A) 

21 M l /n ll 


f (C) a /, 

h 31 A 3 /h 13 


, (c) , . 

n 22 h ll h 33 " h 31 h 13 


h 32 " “ (fl ll h 32 ~ h 12 h 31 ) 


(C) 

23 


li 


(C) 

33 


■ (h ll h 23 ~ ,1 21 ll 13 ) 


h ll h 22 " h 21 h 12 


V 


H -1 = H tC) / |h 


I 




f 


■P'S^ 

,y 







SECTION IV 

MATHEMATICAL FORMULATION 


A. SYSTEM GEOMETRY 
■— > — > — > 

Let i, j, I< be the unit base vectors at the origin of a cartesian reference frame 
XYZ and let be an intersecting pair of baselines in this reference frame, viz 

the AZUSA baseline system. Let R. extend from the origin to either end of B. 

— j» — > J* / __j, I— 

and let R^ R^ extend similarly to either end of Then, from figure 1, 


B L," R 2- R 1, B M = R 4- R 3 


"T* — P 

Then B^ B^ have direction cosines, with respect to X, Y, Z, given by: 


( 1 ) 


cos (B. X) 

L , 

= B L -r/ 

^ 

b l 

• 

i 

cos(b m, X) 

—4 — 

~ b m - i/ 

b m 

■ 

i 

— ► — 

■4 -4 

— > 


— ► 

— ► — ► 


-4 



cos (B. Y) 
*- / 

= B L .j/ 

b l 

• 

j 

c°s (B m Y) 


b m 

• 

j 

• — b “4 

— > —4 

-4 

6 

-V 

— *• -4 

— > — W 

— *> 


— > 

cos (B. Z) 
L , 

= b l v 

b l 

•Ik 

1 

cos (B^j Z) 

= V k/ 

b m 

• 

!< 


( 2 ) 


- r (x) 7 , R (y)“f , R (z)t - r(x)*? , n (y)“? , ptz)"? 
B L ~ B L i+B L j+B L k ' B M- B M 1 +B M j + B M k 


. tion R 7_ n (x)7. , n (y)?r, c (z)7 7_. Q (x) 
then , B^* i - B^ i-i + Bj' + k-i-B^ , 


(v) 

b m ’ ' = b m 


( 3 ) 


— t- — > 


B L’j = B( L y) 


► — > 

R • I ~ R 

°m j d m 


(y) 


— V — > 

B L' k = B L 


-+ -> (z) 

V k=B M 


4-1 


Since i = j = k = 1 , it follows from equations (2) and (3) that 


cos (B, X) 

L / 


cos (B.. X) — 

M, B. 


B ( y) 

■ D ■ 

COS (B, Y) = ,J=- 
L, B. 


cos (B. . Y) — |=f“ 
M, B. 


cos (B, Z) 

L, f 


COS (B, n Z) - =■ 

Mz B, 


X 1 T+ V + V 


R 2 = X ?J + Y 2 j+Z 2 k 


X,i + Yoj + Z„k 


R 4 = X 4 i+Y 4 j+Z 4 l< 


then by equation (1), 


> ^ X Q ” X -J 

cos (B. X) = 

L, B l 


— - 


cos (B^ X) 


X 4 - X 3 


—S' —t> Y ry “ Y .J 

cos (B. Y) = 

L ' b l 


7 “7 

»~-£* £— . A 4 « 

cos (B, Z) = — — 

L ' B L 


cos (B m Y) 


cos(B M, Z) 


y 4 -y 3 


z - z 


4-2 


where: 


b l = b l 


'(X 2 - X ) 2 + (Y 2 - Y,) 2 I- (Z 2 - Z^ 2 



b m = b m = V(X 4 -X 3 ) 2 +(Y 4 -Y3) 2 H-(Z 4 -Z 3 ) 2 


— p - — y , 

Thus, the orientations of with respect to X, Y, Z have been obtained. 


— ■> — ► 


Now let R^ R^ extend from the origin to the centers of respectively. 

Then, from figure 2 , 


B 


L 


B 


M 


R L 2 + R l, R M 2 + R 3 


( 6 ) 


But, substituting from equation (1) 


R 2 " R 1 R 2 + R 1 


R L 2 


+ R 1 - 2 


(7) 


— y 

R P 

*3 


R 4 +R 3 


r m 


2 + R 3 


Let R be the slant range from the origin to the vehicle. 

— >■ — > — j * — ► 

Let L, M be the slant ranges from the vehicle to the centers of B^ B^ respectively 
From figure 3, 


— > — v — ► — > — > — > 

l = r ~ r l, m = r - r m • 


( 8 ) 


Following the reasoning of equations (2), (3), and (4), the direction cosines of 


— — > 

L, M with respect to X, Y, Z are given by: 


-- V 


— > — > i L x 

cos (L, X) = L * = ~T~ 

L • i L 


(9) 




cos (L, Y) 


Ly 


L 
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i 







cos (L/Z) 

**► 

- L / L 

■/! 


-F> V 

cos (M ,X) 

“ M / M 

X 

cos (M ,Y) 

- m /I'm 

y 1 



cos (M ,Z) 

= M / M 
z 


( 


Where L = L i + L j + L k, M = M i H M j + M k 
x y z x y z 


" R L = X ,. i+V L 1 + 2 L k - R M = X M i 1 V +Z M k 


R = Xi + Yj + Zk, then by equation (8) 
cos (L,X) 

— y 

cos (L,Y) 

— — > 
cos (L,Z) 


(X - x l )/l 
(Y - Y l )/L 


cos (M,X) = (X-X m )/M 
cos (M,Y) = (Y-Y m )/M 


(Z - Z. )/L 


cos (M ,Z) = (Z-Z m )/M 


where: L 


L 



(X -X L ) 2 +(Y - Y,) 2 +(Z -Z L ) 2 


( 10 ) 


M = 


M 



- V(X - X L ) 2 + (Y - Y l ) 2 + (Z - Z L ) 2 


Having thus obtained the orientations of L, M with respect to XYZ, we can now 

r "- , K — *► • — ^ — -) »• 

derive the orientations of L, M with respect to respectively: 


— > — > *— >■ 


cos (L ,B^) - L, * B |^/ 


L 


B 


L 


- 4 * — •> — 4 * 


cos (M /Bjyj) = M • B^/ 


M 


B 


M 


( 11 ) 


But 


L'R = L b{ x) I L [i{ y) 'I' L b| z) 
L x L y L l L 


“> = ■" Vm “z B M 

and hence from equations (4) and (9) 


cos (L,B. ) ~ cos (L,X)*cos (B, X) H- cos (L,Y)* cos (B, Y) 

L. L / * L / 


-f cos (L,Z) * cos (B, Z) 

L / 




cos (M,B m ) = cos (M,X)-cos (B m X) + cos (M ,Y) • cos (B^ Y) 
+ cos (M ,Z) ' cos ( ,Z) 


The symbolic notation for these direction cosines is 


A £ cos (L,B^ ) 


m e cos (M , Bjyj) 


The AZUSA system normally measures j£,m,i ,m and a slant range from some 
point (other than the centers of B^ B^) on the system, and an associated range rate, 
To derive this range, let R r extend from the origin to range measuring point and 

— V J 

let R^ extend from the vehicle to the range measuring point. 

I 

Then by the reasoning of equation (8), 


— > -+• 


R R = R - R 5 



,*.s»ciu R is as in equation (8). Upon iho assimmiim. {lint 


R r "XJ-i-Y P MZ r k 
5 D 5 5 


— ► 


R r s (X - X 5 )i *1- (V - Y r )j I (Z - Z 5 )k 


( 14 ) 


let 


To simplify the derivation of the expressions for the rates mentioned previously. 


a L = 

cos (L,X) 

= (X - X.)/L 

a M = 

cos (M ,X) 

= (X - x m )/m 

b L“ 

cos (L,Y) 

= (Y - Y l )/L 

b M 5 

cos (M ,Y) 

“•V — t* 

= (Y - Y m )/M 

c L E 

cos (L,Z) 

= (Z - z l )/l 

C M 5 

cos (M ,Z) 

= (Z - z m )/m 


d^ = cos (B^ X) 

—f nh 

e^ = cos (B^ Y) 


cos (B, Z) 
t - 1 


d M = cos (B M P 

2 m 5 cos (B m P 

fc« - cos (B. > Z) 
M M , 


Then, by equation (12) 


* = a lA + b L e L + 0 L f L 


m a M d M +b M e M +c M f M 


(15) 


Now d^ d^ e^ are all constant and hence have zero derivatives 

with respect to any variable. Consequently, 


= d L a L +e L b L +f L C L 


m d M a M + e M b M + f M c M 


(16) 


» i * • • • 


All that remains then is to evaluate a, b. c. a. fl b.. c., 

L, L, L, M, M, M 


4-T 




w 


Now 




and il is easily verified that 


Sa L 

(1_a L ) 

aa L 

_ " a L b L 

aa L 

_ " a L C L 

SX 

L. 

ay 

L 

a z 

L 

3b L 

_ " a L b L 

ah, 

(1-b?) 

ab L 

_ _b L C L 

ax 

L 

BY 

i 

L 

az 

L 

ac L 

-3 C 

o-,— ^ 

ac L 

_ " b L c L 

a°L 

i 

i-J 

! 

O 

, 

ax 

L 

av 

L 

'az 

L 


♦ 


( 18 ) 




Substitution of equations (17) and (18) into equation (16) yields the desired 
expressions for &, m. 


0 




The range rate R p is the rate of change of the magnitude of the range vector 


in the direction of the velocity vector V = Xi H- Yj + Zk 


T” 

This is so since the vector representing the rate of change of has a maximum 


component in the direction tangent to the flight path at (X,Y,Z), i.e. , in the direction 
•+ — *■ . 

of V. The space rate of change of RjJ is giverfby 


R 


R 


R 


R 


aX 


+ 


a 

-V 

Rd 

| — *• 


R 

j 


R 


R 


3Y 


az 


(19) 


Then the range rate is given by 


R r = V 


R 


R 


V 


( 20 ) 


* 


4* 


10 



i 


which upon expansion is 


SR 


R = 
k R 


R 


SX 


X 


sr 

Ty 


r 


Y + 


B R 


R 


oZ 


( 21 ) 


where R 


R 


R 


R 


Since 


SR 


R __ 


sx 


x-x 


R 


R 


SR 


R 


a Y 


Y~Y 


R 


R 


SR 


R 


sz 


z-z 


R 


R 


Then equation (21) becomes 


(X-XJX 


(Y-YJY (Z-ZJZ 
n - 'j . 5 , 5 

R R R n Rn H R ‘ 


l R 


R 


R 


NOTE: The derivation of die equations used to transform geographic data 

to cartesian data can be found in Reference 1. 


B. MATHEMATICS OF THE COOP DETERMINATION 


GDOP determination is the problem of dolomnnincj to what degree errors in param- 
eter measurements are carried through non-orthogonal transformations. The solution 
to the problem treated here is very specialized to simplify computations. In particular, 
it will be assumed that velocity contributes no error to position, the effect of which 
assumptions will be seen further on. The general development of a GDOP determination 
scheme will not be presented hero. Also, see Reference 2. 


Let $., be the position covariance matrix, which is a diagonal matrix consisting 
1 2 2 2 

of the variances a o <j p (this presumes independence of errors for z , m, R). 

Let B be the matrix of partial derivatives of z , m, R with respect to X,Y,Z. Then tl e 
position covariance matrix S 2 for X,Y,Z (i.e. , the matrix indicating the extent of prop- 
agated error) is given by 



B T S” 1 B 


-1 


(23) 


/ 

t 


Then the parameter variances for X,Y,Z are the square roots of the diagonal 
elements of S£ / i .e . , 


a 


x 


.( 2 ) 

ll 


a 


= S 


( 2 ) 

n 


a 


= S 


( 2 ) 

33 


(24) 


c\ 

w 


where 


s (2) 

IJ 


3 

ij 


The case for velocity is similar, except that the partial derivative matrix B 

• • • « • 1 

contains the derivatives of % , in , R by X, Y, Z, and in the place of S 1 the 

covariance matrix, S- for z, m, R is substituted, yielding the covariance matrix, 
• • • ^ ^ 

S 4 for X, Y, Z. 


4-12 


o 


Now from equations (15) and (18) in paragraph A, 


L d 


L '!■ 


L e 


L 4- 


'L f, 


( 25 ) 


(1-a? ) , 

L d, 


a L b L e L 


a L°L f L 


" a L b L d L , (1-b L )c L 

" L " L 


b L c L f L 


_a L c lA b L°L f L (1 ' c L )f L = C, 
l + L ‘ L 


(1 “ a M )d M 


a M b M e M 


a M c M f M r A, 


‘ a M b M d M j. (1-b M )e M 
M M 


h p f 

M M 


9 

c d hce O -c )f 

”M M M M M M , u M M 

M ” M M 


and also from the statement preceding equation (22) in paragraph A 


X-X 


3R R _ Y - Y 5 


4-13 


The computed values of equations (25) mid (26) are substituted into equation 
(23) giving , 


( 27 ) 


S 2 = 


a l a m a r 


b l B Wi b r 

C L C M C R 


"2 n 

ff A • 0 


0 


0 0 ° r 


a l b l c l 


a m b m c m ’ -1 
a r b r c r 


When the matrix product in the brackets is expanded, a matrix A is obtained 
whose elements are 


A R 
a R 


(28) 


A L B L 


A M B M 


a r b r 


a l c l 


a m c m 


a r c r 
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0 


A 21 

A 12 

A 23 

_ b l c l 
•? 

r — i 
< 

A 13 

A 32 

A 23 

The inverse of A, 


, b m c m . b r c r 

H* * — "r — ■=■ — 


a 


m 


a 


R 


method. 

In general it is not true that 
* 

d H d 


P 
dX 


p 


ax 


for any parameter p, but in the present situation it is true that 

SR, 


dil 5 A 


Bm 


dm 


R 


dR 


R 


dX 


dX 


dX dX 


X 


dX 


, etc. 


To see that this is so, recall that (from equation (16) paragraph A) 


d L a L + e L b L + f L c L 


m d M ? ‘IVl + e M b M + f M c M 


and that from equations (17) and (18) this yields 


4-15 


i 

i 

i 

i 


a -"L M t. 


\ b L e L VlA y 

L~ ~L 


(29) 


“ a L b L d L 


(i-b^L 


b L c L f L. x, 
L 


-a L c L cl L b L c L°L , (1 ' c L )f L ; 

L L L L 


(1 ' a IVl )d M 


a M b M e M 


a M c M f M y 
M X 


' a M b M d M 


Q-b M )e M 


b M C M F M ' 
M 


j " a M C M d M b M c M e M . <1-c M )f M ; 
M " M M 


Upon differentiating equation (29) by X, Y, Z, it is seen that 


Q - a L )d L 


a L b L e L a L C L f L 


(30) 


" a L b lA , ( 1 ' b L? b L C L f L 

— — ~r 


-a L c L d L 


b, c, e, 


d-b 2 ) 
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dm 

BX 


JVl M 

M 


a M b M°M 

M 


a M fM f M 
M 


Bm 

BY 


' a M h M d M 

M 




M 


[) (• f 

VjA 

M 


B m __ 


‘ a M C M d M 


b M c M e M 


(1 “ C M )f M 


BZ 


M 


M 


M 


Upon examination of equations (30) and (24), it is seen that 


= A 


BX 


M 


ju 

BY 


B, 


B.e _ 


BZ 


(31) 


Bm __ 


A 


BX 


M 


im = B 


BY 


M 


Bm 


r = C 


BZ 


M 


Recalling equations (22) and (26), 


BR 


R 


(X-X,) 

5 _ 


BX 


A 


R 


R 


R 


Bf^ 

BY 


(Y-YJ 

5 


= B 


R 


R 


R 


(32) 


BR 


R 


(Z-Zr) 


B Z 


R 


'R 


R 


Consequently, for computation purposes, once B has been computed for position, 
it is saved and used in the velocity solution, since B is the same in both cases. 


J 
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/ 




' f wvn** w v 


'Mr* » i ; 


A computation is also made which propagates X,Y,Z errors into the magnitude 
V of the velocity vector V « Xi H* Yj T 7-k. Note that this is a just-dotemiined 


* • • 


condition since X, Y, Z completely determine V. In such a case the covariance matrix 
S 5 for V is given by 


C s 4 C 


(33) 


where S* is the covariance matrix of X, Y, Z, as in tlie remarks preceding equation 

* r f I » 

(24), and C is the matrix of partial derivatives of V ~ V with respect tc X, Y, Z, 


It will be seen that 


and thus C = r. 


X Y Z 
V V V 


dV = Y 
BY V 


. Then equation (33) becomes 


(34) 


X Y Z 
V V V 


ry . ry • • O’ * * 

J x u xy >:z 


2 

o • • o • 
xv v 


a • • 

yz 


( 35 ) 


a . . a . . 
xz yz 


which on expansion becomes 


(36) 


S, =[Vl [x 2 ^ + S<2 ° ?: + + 2XYa • • + 2XZa- • + 2YZa- •] 

5 L,2 L x y z xy xz yz J 


4-18 



Since has only one element, 


2.2 , •, 2,2 , ; 2_2 


V 


X a. + Y a. H* Z fc o. +2XYa- + 2XZo*. -1* 2YZ a -- 

_x y z xy xz ^az 

V 


(37) 
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SECTION V 

COMPUTER PROGRAM OPERATING INSTRUCTIONS 


A. PROGRAM DESCRIPTION 


1 . 

Program Identification 

AZGD2 

2. 

Computer 

GE-635 

3. 

Program Library Number 

1016 

4. 

Type of Coding 

Fortran IV 

5. 

Core Storage 

6666 


B. USAGE AND LIMITATIONS 

This program computes GDOPS (accuracies) for a standard KSC theoretical 
trajectory based on an AZUSA (or AZUSA-like) system; that is, a tracking system 
with two intersecting baselines. 


5 * 


PROGRAM INPUT FORMAT 

t 

1. Tape Input 

Standard KSC trajectory tape, or with the following format will be used: 
Word No , Entry 

1 Time (sec) 

2-4 X, Y, Z (mtrs) 

X, Y, Z (mtrs/sec) 


5-7 


Card Input 
a* Deck Setup 


SNUMB 


$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 


IDENT 

OPTION FORTRAN 


PROGRAM DECK 


EXECUTE 


LIMITS 
TAPE 03 
FFILE 03 
TAPE 04 
INCODE IBMF 
DATACARDS 
ENDJOB 


Operations Supplied (OS 
Requester Supplied (RS) 
. S 
OS 


RS 

RS 

RS 

RS 

RS 

RS 

RS 

RS 


***EOF 


b. Card Formal; 


Card 

Columns 

In forma lion 

Format 

Remarks 

1 

1 


$* 

Control card 


8-12 


IDENT 



16-20 

Program ID 

AZGD2 



21 


/ 

Comma 


22-27 

Requester's 
control number 

xxxxxx 



28 


t 

Comma 


29-31 

Requester's 

initials 

XX 



32 


l 

Comma 


33-34 

Project code 

XX 



35 





36-37 

Rerun code 

XX 



38 





39-48 

Work order 
number 

xxxxxxxxxx 


2 

1 


$ 

Control card 


8-14 


EXECUTE 


3 

1 


$ 

Control card 

• 

8-13 


LIMITS 



16- 

Maximum run 
time / hr/100 




Comma 


5 '"4 


Card 

Columns 

In formal ion 

Format 

3 


Core storage 

8000 

(Cont) 


Maximum 
number of lines 
to print 

/ ; 

4 

1 


$ 


8-11 


TAPE 


16-17 

Kile code 

03 


18 


; 


19-21 

Logical unit 
designator 

AIR 


22-23 

i 

! 1 


24-28 

File serial 
number 

i 


29-30 


/ / 


31-42 

File name for 
external use 


5 

1 


$ 


8-12 

* 

FFILE 


16-17 

File code 

03 


f 


Remarks 

Commas 

*10000 

Control card 

Zero-three 

Comma 

If input tape is to 
be returned to 
requester, replace 
R with D. 

Commas 

99999 if unknown 
Commas 

Information on tape 
label 

Control card (only 
present for tape 
input) 

Zero-three 


18 


Comma 


Card 

Columns 

Information 

Format 

Remarks 

5 

19-24 


NSTDBL 

No standard labels 

(Cont) 

25 


/ 

Comma 


26-34 


BUFSiZ/54 

54 words per record 

* 

35 


/ 

Comma 


36-44 


FIXLNG/54 

54 words per record 


45 


! 

Comma 


46-51 


N0SRLS 

No block serials 

6 

1 


$ 

Control card 


8-11 


TAPE 



16-17 

File code 

04 

Zero -four 


18 


/ 

Comma 


19-21 

Logical unit 

B1R 

If output tape is to 



designator 


be saved, use S or 





D for R . 

If output tape is needed the following fields must be used on 

Card 5: 

o 

2 2 -2d 


/ / / / 

Commas 


26-37 

File for 
external use 

t 


7 

1 


$ 

Control card 


8-13 


INC0DE 



16-19 


I BM F 


8 

1-5 

Program ID 

AZGD2 

Comment card 


6-80 

Anything 




Card 

9 


10 


11-16 


Columns Ini urination 


Format 


Roma 


1 Coordinate X 
system of 
output 


0~oulpuL in coordinate 
system 

/ X-'downrange \ 

Y=u|) 


\ Z-crossrange / 
1-output in Apollo 
Saturn Standard 
Coordinate System 
I Z-downrange 
X=up 

\ Y-crossrange 


1-15 

a of spheroid ±0 .XXXXXXXXE±XX 

Meters 

16-30 

b of spheroid ±0 .XXXXXXXXE±XX 

♦ 

Meters 

31-45 

Azimuth of ±0 .XXXXXXXXE±XX 

X-axis 

Degrees 

1-9 

. 


10-12 

cp j z Degree ±XX 

Trajectory origin 


j-0: trajectory 
origin 

Latitude 

* 

j“l:Ll 

Coordinate 

L Baseline 


j=2:L2 

Coordinate 

Latitude 


j=3:Ml 

Coordinate 

M Baseline 


j=4:M2 

Coordinate 

Latitudes 


jj—5: Origin 
the R parameter 

R origin 
latitude 


5-7 


Card 

Coi limns 

Information 

11-16 
(Co nt) 

14-15 

T' j , Minute 


17-24 

v j , Second 


27-30 

Kit Degree 


32-33 

A j, Minute 


35-42 

A j, Second 


45-56 

hj, Meters 


57 



58 

‘Spheroid 


59 



60-80 

ID of loca- 
tion j 

i 


Example: 

(Card 11) 
(Card 12) 
(Card 13) 
(Card 14) 
(Card 15) 
(Card 16) 


17 1-15 

16-30 a m 

31-45 ct r 
46-60 cj£ 
61-75 


Formal , 

XX 

xx.xxxxx 

d'-XXX 

XX 

xx.xxxxx 

iXXXX.XXXXXX 

X 


Pad Coordinate 
LI Coordinate 
L2 Coordinate 
Ml Coordinate 
M2 Coordinate 
R Coordinate 

cards and must be in the order 

±O.XXXXXXXXE±XX 

±0.XXXXXXXXE±XX 

±O.XXXXXXXXEiXX 

±0.XXXXXXXXE±XX 

±0.XXXXXXXXE±XX 


Cards 11-16 are ail coordinate 


Remarks 

Longitude 

Longitude 

* 

Height 

F, C, W or K 


shown above. 

L-direction cosine 
accuracy 

M-direction cosine 
accuracy 

Range accuracy (mtrs) 
L-rate accuracy 
M-rate accuracy 


Card 

Columns 

Information 

Format 

Remarks 

18 

1-15 

°R 

0 .XXXXXXXXD XX 

Range-rate accuracy 
(mtrs/sec) 

19 

1-15 

Start time 

*O.XXXXXXXXE*XX 

Second 


*16-30 

Stop time 

±O.XXXXXXXXE:tXX 

Second 

20 

1 


$ ' 

Control card 


8-13 


ENDJOB 


21 

1-6 


***EOF 

End of file 


NOTE: The program normally calls for a Fortran binary output tape on logical 
File 04. If the tape is not needed, leave Columns 22-80 blank on 
Card 6. 


PROGRAM OUTPUT FORMAT 


The output data are a printer output listing. 



E. SAMPLE TEST CASE 
1. Input 


a. Test One 
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s. 


V, 


COORDINATE SYSTEM 

)<oDoHNR 4NG£ (DR}; Y*UPiZ»eROSSRANOE (90>DR! 


2. Output 

a. Tape Output Format 


Word No * 
1 

2-8 


Entry 

Time (sec) 

cr v o,. o 7 (mtrs) 
x , y t 4 

a • a- a- o- (nitrs/sec) 
x, y, z, v 


b. Test One 


/ 


16 


TIME 

sec 

0, 


S* ' vv 

M/S 

O.32O71069B 

04 

SIG 

z 

M 

0,969846516 

02 

L 5 DZ 

M/S 

0 . 24 ?4935 OE 

01 

S l G 

X 

M 

0,803024386 

05 

SIG DX 

M/S » 

0.220706286 

04 

* SIG 

V 

M 

0.2B355843E 

03 

SIGYDZ 

M/S 

0,706952816 

01 

TIME 

see 

o.idooooqoe 

01 

SIG VV 

M/S 

0,221060076 

04 

SIG 

i 

M 

0 * 1 0199760E 

03 

SIG DZ 

M/S 

0 , 2530901 66 

01 

SIG 

X 

M 

0,883643226 

05 

SIG DX 

M/S, 

0 , 221081056 

04 

SIG 

V 

M 

0,2907 1130E 

03 

SIGYDZ 

M/S 

0 , 72730152E 

01 

* AM 

time 

SEC 

0.20000000E 

01 

SIG VV 

M/S 

0,221140576 

04 

SIG 

l 

M 

0 , 11739044E 

03 

SIG DZ 

M/S 

0,293522136 

01 

SIG 

X 

M 

0.B8451439 E 

05 

SIG DX 

M/S. 

0,221148996 

04 

SIG 

Y 

M 

0, 31274728E 

03 

SIGYDZ 

M/S 

0,781903705 

01 

‘TIME 

SEC 

0.30000000E 

01 

SIG VV 

M/S 

0.221676316 

04 ' 

SIG 

Z 

M 

0.14429137E 

03 

SIG DZ 

M/S 

C .360902786 

01 

SIG 

X 

M 

0 , B0627494E 

05 

SIG DX 

M/S, 

0.221676326 

04 

SIG 

Y 

M 

0.35U6767E 

03 

SIGYDZ 

M/S 

0.87831 2606 

01 

time 

SEC 

0,40000000E 

01 

SIG VV 

M/S 

0 . 22257035E 

04 

“SIG 

Z 

M 

D.1B355324E 

03 

SIG DZ 

M/S 

O.45079864E 

Oi 

SIG 

X 

M 

0 /89046787E 

05 

SIG DX 

M/S, 

0,222570236 

04 

SIG 

Y 

M 

0, 40753543E 

03 

SIGYDZ 

M/S 

0.101859646 

02 

time 

SEC 

0 .500000Q0E 

01 

9 I G VV 

M/S 

0 . 22482259E 

04 

SIG 

1 

M 

0 , 23452763E 

03 

SIG DZ 

M/S 

0.591238516 

01 

SIG 

‘X 

M ' 

"0 .89180629E 

05 “ 

S I G DX 

M/S, 

0,224022406 

04 ~ 

SIG 

Y 

M 

0 , 40OO2524E 

03 

SIGYDZ 

M/S 

0.12101078E 

02 

TIME 

SEC 

0 , 6 0 0 000 0 0E 

01 

SIG VV 

M/S 

0,224765876 

04 

SIG 

Z 

M 

0.3044B053E 

03 

5 I G DZ 

M/S 

0,749948196 

01 

SIG 

X 

M 

O.91255270E 

05 

SIG DX 

M/S, 

0.22476564E 

04 

"'SIG 

V 

M 

' O.50264159E 

03 

SIGYDZ 

M/S 

0.14350528E 

02"“ 

TIME 

SEC 

0,700000006 

01 

SIG VV 

M/S 

0,225802195 

04 

SIG 

1 

M 

0 , 3803532 2E 

03 

SIG DZ 

M/S 

0,945127446 

01 

SIG 

X 

M 

0 , 90870768E 

05 

SIG DX 

M/S, 

0 , 225801966 

04 

SIG 

Y 

M 

0.6895073 3 E 

03 

SIGYDZ 

M/S 

0 , 17133178E 

02 

time 

SEC 

0 , BOOOOOOOE 

01 

SIG VV 

M/S 

0 , 2374A598E 

04 

SIG 

Z 

M 

0 . 486207926' 

03 

SIG DZ 

M/S 

0 ,123043936 

02 

SIG 

X 

M 

0 ,930424720 

05 

SIG DX 

M/S, 

0 . 2374fl5flOE 

04 

SIG 

Y 

M 

0 . B4 459590E 

03 

SIGYDZ 

M/S 

0.21 37391 7E 

02 

‘TIME 

SEC 

0 .9QOOOOOOE 

Cl 

SIG VV 

M/S 

0.243176406 

04 

SIG 

7 

A* 

M 

0 , 6076Q546E 

03 

SIG DZ 

M/S 

0 . 153834706 

02 

SIQ 

X 

M 

0.960482326 

05 

SIO DX 

M/S, 

0 , 2 *’3176356 

04 

SIG 

Y 

M 

0 , 102049105 

04 

SIGYDZ 

M/S 

O.25836010E 

02 

T l ME 

SEC 

0.10000Q00E 

02 

SIG VV 

M/S 

0 , 255252178 

04 

' SIG 

1 

M 

0 ,787610695 

03 

SIG DZ 

M/S 

0 ,194600286 

02 

StG 

X 

h 

0 .103266495 

06 

SIG DX 

M/S, 

0.253252376 

04 

SIG 

Y 

M 

0,128755035 

0 4 

SIGYDZ 

M/S 

0,310233596 

02 


i 

i 

i 


c. Test Two 


T I ME 

SEC 

0, 


SIG 

vv 

M/S 

0.22071069E 

04 

SIG 

X 

M 

0 * 969B4631S 

02 

SIG 

DX 

M/S 

0.24249350E 

01 

SIG 

V 

M 

0,803024388 

03 

SIG 

DY 

M/S, 

0 i 22O7062BE 

04 

SIG 

z 

M 

0,283550438 

03 

3 I G 

DZ 

M/S 

0,700992011: 

01 

TIME 

SEC 

0 , 100000008 

01 

SIG 

VV 

M/S 

0.22108007E 

04 

SIG 

X 

M 

0 ,10193^606 

03 

SIG 

DX 

M/S 

0.29909016E 

01 

SIG 

V 

M 

0 .08364322E 

03 

SIG 

DY 

M/S , 

0 , 2210H185E 

04 

SIG 

z 

M 

0,290711306 

03 

SIG 

DZ 

M/S 

0 , 72730152E 

01 

TIME 

SEC 

0,200000006 

01 

SIG 

VV 

M/3 

0 , 22114857E 

04 

SIG 

X 

M 

0,117398446 

03 

SIG 

DX 

M/S 

0 .29352213E 

01 

SIG 

Y 

M 

0,804314396 

OS 

SIG 

DY 

M/8, 

O.22114099E 

04 

SIG 

Z 

M 

0,31274/206 

03 

SIG 

DZ 

M/S 

0 , 78190370E 

01 

TIME 

SEC 

0,300000006 

01 

9 I G 

VV 

M/S 

0 .22167631E 

04 

SIG 

X 

M 

0 .14429137E 

03 

SIG 

DX 

M/S 

0 , 36O9O270E 

01 

SIG 

Y 

M 

0 , 886274946 

05 

SIG 

DY 

M/$, 

0.22167632E 

04 

SIG 

Z 

M 

0,351167671 

03 

8IG 

DZ 

M/S 

0.87831268E 

01 

TIME 

SEC 

0.400000006 

01 

SIG 

VV 

M/S 

0.22297035E 

04 

'SIG 

X 

M 

0 , 183558246 

03 

SIG 

DX 

M/S 

0 , 43079B64E 

01 

SIG 

Y 

M 

0 , 890467876 

03 

SIG 

DY 

M/S, 

0.22257023E 

04 

SIG 

Z 

M 

0 , 40753543E 

03 

SIG 

DZ 

M/S 

0.10185964E 

02 

TIME 

SEC 

0,500000006 

01 

SIG 

VV 

M/S 

0 , 22482259E 

04 

SIG 

X 

M 

0,234527636 

03 

SIG 

DX 

M/S 

0 .39123891E 

03 

‘ S I Q 

Y 

M 

0.89180629E 

03 

SIG 

DY 

M/S, 

O.2240224OE 

04 

SIG 

Z 

M 

0,480025246 

03 

SIG 

DZ 

M/S 

0 .12101078E 

02 

TIME 

SEC 

0 , 60 0000006 

01 

SIG 

VV 

M/S 

0.22476567E 

04 

SIG 

X 

M 

0,304480536 

03 

SIG 

DX 

M/S 

0 , 74994S19E 

01 

SIG 

Y 

M 

0 , 91255278E 

03 

SIG 

DY 

M/S, 

0 , P2476564E 

0 4 

SIG 

Z 

H 

0,582641396 

0 3 

SIG 

DZ 

M/ 8 

0 , 1435052BE 

02 

TIME 

SEC 

0 , 70000Q00E 

01 

SIG 

VV 

M/S 

0.225B0219E 

04 

SIG 

X 

M 

0,380353126 

03 

5 I G 

DX 

M/S 

0 , 94512744E 

01 

SIG 

Y 

M 

0 .9087076SE 

05 

SIG 

DY 

M/S, 

0.22580196E 

04 

SIG 

1 

M 

0 , 68958733E 

03 

SIG 

DZ 

M/S 

O.1713S170E 

02 

TIME 

SEC 

0,800000008 

01 

SIG 

vv 

M/8 

0 , 2374O09OE 

04 

SIG 

X 

M 

0 ,486207926 

03 

SIG 

DX 

M/S 

0 , 12304393E 

02 

SIG 

Y 

M 

0 ,930424728 

03 

SIG 

DY 

M/S, 

0 , 237 4B380E 

04 

SIG 

1 

M 

0,844595906 

03 

BIG 

DZ 

M/S 

0 , 21373917E 

02 

'TIME 

'SEC 

0 ,9 0 0 000006 

01 

SIG 

vv 

M/S 

0.24317640E 

04 

SIG 

X 

M 

0,607685466 

03 

SIG 

DX 

M/S 

0,13389470E 

02 

'SIG 

Y 

M 

0,960482326 

05 

SIG 

DY 

M/S , 

0,24317639': 

04 

SIG 

1 

M 

0,102049108 

04 

SIG 

DZ 

M/S 

0 , 25O36018E 

02 

TIME 

SEC 

0,100000008 

02 

8 I G 

vv 

M/3 

0 .29529217E 

04 

SIG 

X 

M 

0,787610696 

03 

SIG 

DX 

M/S 

0 ,194630286 

62 

3 1 G 

Y 

M 

0,103266496 

06 

SIG 

DY 

M/9, 

0 , 2592S237E 

04 

SIG 

Z 

M 

0,128755038 

0 4 

SIG 

DZ 

M/8 

0,3182935 9 E 

02 


5-17/18 


SECTION VI 

COMPUTER PROGRAM FLOWCHART 


START 


SET LINE s 13 
(LIME CONTROL) 



READ ID 
COMMENTS 


READ ORIGIN 

\ , U , in DEG 
ooo 

MIN. SEC. MTRS. 


PRINT GCOGR 
TITLES; PRINT 
,ORG!N 
LOCATION 


IS 

ID AZGD2 
? 


PRINT MESSAGE, 
PROGRAM ID 
(AZGD2) DOES 
NOT CHECK 
WITH CARD 
ID (XXXX). 



/ 

READ 

AND 

> o, b, 

AZ. 

l 

7 

f 


PAGE 6-2 


CALL DEGRAD 
TO CONVERT 
<!> 0 > A 0 , TO 


RADIANS 



DO 

J ~ 1, 5 


HALT 

(EXIT) 


COMPUTE e 2 


SLEW ONE 
PAGE, PRINT 
ID, COMMENTS; 
PRINT a, b, AND 
AZ 



READ BASELINE 
END POINT 6 j,Aj 
H J( IN DEG. MIN. 
SEC. MTRS. 


PRINT GEOGR 
TITLES; PRINT 
LOCATION OE 
POINT J 




CONVERT AZ 
TO RADIANS 


PAGE 6-2 






PACK 6-1 



















DEGRAD ROUTINE 


GEOXYZ ROUTINI 



NOTE! 

A ~ DEGREES 
B « MINUTES 
C «, SECONDS 
D s RADIANS 



6-4 








PHIGAM ROUTINL' 















VARCOV ROUTINE 


MINVRT ROUTINE 
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